Deconfinement of neutron star matter within the Nambu-Jona-Lasinio model 
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We study the deconfinement transition of hadronic matter into quark matter under neutron 
star conditions assuming color and flavor conservation during the transition. We use a two-phase 
description. For the hadronic phase we use different parameterizations of a non-linear Walecka model 
which includes the whole baryon octet. For the quark matter phase we use an SU (3) / Nambu-Jona- 
Lasinio effective model including color superconductivity. Deconfinement is considered to be a first 
order phase transition that conserves color and flavor. It gives a short-lived transitory colorless- 
quark-phase that is not in /3-equilibrium, and decays to a stable configuration inr~ T wea k ~ 10 _8 s. 
However, in spite of being very short lived, the transition to this intermediate phase determines 
the onset of the transition inside neutron stars. We find the transition free-energy density for 
ff") ' temperatures typical of neutron star interiors. We also find the critical mass above which compact 

stars should contain a quark core and below which they are safe with respect to a sudden transition 
to quark matter. Rather independently on the stiffness of the hadronic equation of state (EOS) we 
find that the critical mass of hadronic stars (without trapped neutrinos) is in the range of ~ f .5 - 
1.8 solar masses. This is in coincidence with previous results obtained within the MIT Bag model. 
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■ I. INTRODUCTION 

The core of compact stars can reach densities that are several times larger than the saturation density of nuclear 
matter. In such extreme conditions the baryons get so compressed that they start to overlap and can produce a 
deconfined gas of quark matter. An important characteristic of the deconfinement transition in neutron stars, is that 
just deconfined quark matter is transitorily out of equilibrium with respect to weak interactions. In fact, depending 
on the temperature, the transition should begin with the quantum or thermal nucleation of a small quark-matter drop 
near the center of the star. On the other hand, the flavor composition of hadronic matter in /3-equilibrium is different 
from that of /3-stable quark-matter drop. Roughly speaking, the direct formation of a /3-stable quark-drop with N 
quarks will need the almost simultaneous conversion of ~ N/3 up and down quarks into strange quarks, a process 

I-"" ' 2N/3 

■ which is strongly suppressed with respect to the formation of a non /3-stable drop by a factor ~ G Fcrmi . For typical 
Q\ , values of the critical-size /3-stable drop (N ~ 100-1000 hi) the suppression factor is actually tiny. Thus, quark flavor 

must be conserved during the deconfinement transition p], 0, [H, 0, [H, H 0] • The main consequence of this condition 
is that the density of the transition is higher than it would be if the direct formation of /3-stable quark-drops were 
possible. This is easy to understand, since the Gibbs free energy per baryon of /3-equilibrated quark matter is always 
smaller of that of the non- /3-equilibrated state [f|. 

Due to the uncertainties in the knowledge of the state of matter at the densities of interest, studies of the decon- 
finement transition are usually based on the extrapolation to higher densities of an hadronic model valid around the 
nuclear saturation density po, and the extrapolation to ~ po of a quark model that is expected to be valid only for 
high densities. Within this kind of analysis the (in general) different functional form of both EOSs, induces the phase 
transition to be first order. Notice that from lattice QCD calculations there are indications that the transition is 
actually first order in the high-density and low-temperature regime, although this calculations involve temperatures 
that are still larger than those in neutron stars, and do not include the effect of color superconductivity In a 
previous work [6j, the deconfinement transition has been analyzed within the frame of the MIT bag model paying 
special attention to the role of color superconductivity. In the present paper we shall analyze the deconfinement 
transition employing the Nambu-Jona-Lasinio model in the description of quark matter. For the hadronic phase we 
shall use a model based on a relativistic Lagrangian of hadrons interacting via the exchange of a, p, and to mesons 
@. For simplicity, the analysis will be made in bulk, i.e. without taking into account the energy cost due to finite 
size effects in creating a drop of deconfined quark matter in the hadronic environment. 

This article is organized as follows. In Sec. [IT] we briefly outline the non-linear Walecka model used to describe 
the hadronic phase. In Sec. [TIT] we provide some details of the Nambu-Jona-Lasinio model we employ to describe the 
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TABLE I: Parameters of the hadronic equation of state. For each parametrization we give the maximum mass M ma x of a 
hadronic star. 

quark matter phase paying special attention to the conditions of color and flavor conservation. In Sec. IIVI we present 
and discuss our numerical results. Finally, in Sec. [V] a summary with some conclusions is given. 



II. THE HADRONIC PHASE 



For the hadronic phase we shall use a non-linear Walecka model (NLWM) 0, [T(| [TTJ which includes the whole 
baryon octet. The Lagrangian of the model is given by 

C = Cb+Cm+Cl, (1) 
where the indices B, M and L refer to baryons, mesons and leptons respectively. For the baryons we have 



Cb = ^ tpB 



7 M (id p - g^B w p - g P B f • Pn) - {m B - g<jB a) 



(2) 



with B extending over the nucleons N = n, p and the following hyperons H = A, £ + , £ , £ ,5 , and S . The 
contribution of the mesons a, ui and p is given by 
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where the coupling constants are 



g a B = XaB g<y, gu)B — XujB 5", ffpS = X p B ffp 



(3) 



(4) 



where x a B , i u b and x p b are equal to 1 for the nucleons and acquire different values for the other baryons depending 
on the parametrization (see below). The leptonic sector is included as a free Fermi gas which does not interact with 
the hadrons, i.e. 



Cl = ^ ipi (ip - mi) (5) 
i 

There are five constants in the model that are determined by the properties of nuclear matter, three that determine 
the nucleon couplings to the scalar, vector and vector-isovector mesons ga/rria, guj/m^, g p /m p , and two that determine 
the scalar self interactions b and c. It is assumed that all hyperons in the octet have the same coupling than the 
A. These couplings are expressed as a ratio to the nucleon couplings mentioned above, that we thus simply denote 
x fj , x^ and x p . In the present work we use three parameterizations for the constants. One of them is the standard 
parameterization GM1 given by Glendenning-Moszkowski |9[, as shown in Table I. This parametrization employs 
"low" values for x a , x^ and x p . Larger values of these couplings make the EOS stiffer and increase the value of the 
maximum mass of hadronic stars to values above than 2 Af Q , see Table I. 



III. THE QUARK PHASE 

In order to study the just deconfined quark matter phase we use an SU (3) / NJL effective model which also includes 
color superconducting quark-quark interactions. The corresponding lagrangian is given by 

8 

C = ip (ip - m) ip + G ^ Uip r a %j)) 2 + (-0 ij 5 r a ip) 2 

a=0 
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+2H ^ [(V> HstaXa' ipc) (V>c H5Ta^A' i>)] (6) 

A.A'=2,5.7 

where rh = diag(m„, rrid, m s ) is the current mass matrix in flavor space. The matrices, and Xi with i = 1, ..,8 are 
the Gell-Mann matrices corresponding to the flavor and color groups respectively, and To = a/2/3 1/. Finally, the 
charge conjugate spinors are defined as follows: ipc = G ip T and ipc = ip T C, where tp = ip^ is the Dirac conjugate 
spinor and C — i~f 2r y°. 

The lagrangian in Eq. © leads to local chirally invariant current-current interactions in the quark-antiquark and 
quark-quark channels. The latter is expected to be responsible for the presence of a color-superconducting phase in 
the region of low temperatures and moderate chemical potentials. Here for simplicity we do not include flavor mixing 
effects. The values of the quark masses and the coupling constant G can be obtained from the meson properties in the 
vacuum. On the other hand, an estimate of H/G can be obtained from Fierz transformation of the one-gluon-exchange 
interactions in which case one gets H/G = 0.75, which is the value we will use hereafter. 

To be able to determine the relevant thermodynamical quantities we have to obtain the grand canonical thermody- 
namical potential fl(T,fif c ) at finite temperature T and chemical potentials fif c . Here, / = (u,d,s) and c = (r,g,b) 
denote flavor and color indices respectively. For this purpose we start from Eq. ([5]) and perform the usual bosonization 
of the theory. This can be done by introducing scalar and pseudoscalar meson fields a a and 7r a respectively, together 
with the bosonic diquark field Aa- In this work we consider the quantities obtained within the mean field approxi- 
mation. Thus, we only keep the non-vanishing vacuum expectation values of these fields and drop the corresponding 
fluctuations. For the meson fields this implies a = a a r a — diag(a„, (Td,a s ) and ir a = 0. Concerning the diquark mean 
field, we will assume that in the density region of interest only the 2SC phase might be relevant. Thus, we adopt the 
ansatz A5 = A7 = 0, A2 = A. Integrating out the quark fields and working in the framework of the Matsubara and 
Nambu-Gorkov formalism we obtain 



^mfa(T, H fc , CTf, |A|) = -- 



d 3 k 

(2^)3 



In det 



_i_ 1 l^ 1 _i_ ~2 , 2\ , l^l 2 



where 



g-i _ ( # - M + 70A -A75T2A2 
I A*75T 2 A 2 #-M-7o£ 



(7) 



(8) 



Here, we have used fc = (\[2n + l)irT\k\, M = diag(M u , M4, M s ) with Mf = mj + 07, and fi = diag({/i/ c }) in 

flavorOcolor space. The determinant of this 72 x 72 matrix can be calculated analytically if M u = Mg. Thus, in what 
follows we will use the approximation m u = md = m and a u — ad = a which implies M u = Md = M. A detailed 
procedure for the calculation of the determinant can be found in Ref.[ll, [IH, 13 • The resulting contribution to the 
thermodynamical potential is 



j 1 r d 3 fc f A d 3 k ■J-y 

■ 22 j (2^)3 ^dct^- 1 (k)/T] = 2 J J^y J2 u&iiVi) + constant 



2 ^ J (2tt) 3 l w/ ' J (2tt) 3 ^ 

n J y 1 J \ } a —i 

where A is the cut-off of the theory and Lu(x,y) is defined by 

u(x,y) = - [x + rin[l + e-( :E -^/ T ]+rin[l + e-^ + ^/ T ] 

with 



xi t2 = E , £3,4,5 = E s 
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A 2 



(9) 



(10) 



2/1 = Mub , 2/2 = Vdb , 2/3 = fJ-sr , 2/4 = Msg , Z/5 = Vsb , 2/6, 7 = 7, > 2/8,9 = 7, 



(11) 



where E = \JP + M 2 , E s = J k 2 + M 2 
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5.5 
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set 2 


5.5 


110.05 


631.4 
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TABLE II: The two sets of NJL parameters. 



Finally, we include the thermodynamical potential for non-interacting ultra-relativistic electrons 



u 4 u 2 T 2 7tt 2 T 4 
O e = — % - ^— - , (12) 

12tt 2 6 180 



where \i e is the electron chemical potential. Therefore, the total potential for the quark matter electron system is 

12 



1 f 1 |Ar 

Cl(T,{n fc },fi e ,a,a s ) = — ^ dkk 2 Y^u{xi,yi) + ^(<rl + <7d + ° 2 s) + -^f-^c + n e - 



(13) 



Here we have subtracted the constant f2 vac in order to have a vanishing pressure at vanishing temperature and 
chemical potentials. From the grand thermodynamic potential Q we can readily obtain the pressure P = — f2, the 
number density of quarks of each flavor and color 

= ~Wo = £ dk k2 Wf c »>) < (14) 

and the number density of electrons 

0/J, e 

The corresponding number densities of each flavor, rif, and of each color, n Cl in the quark phase are given by 



(15) 



n f = J2 n f c ' Uc = n f c - ( 16 ) 

/ 

Finally, the baryon number density ns reads 

n B = \^^n fc = ^(n u + n d + n s ) (17) 

fc 

and the Gibbs free energy per baryon is 

ffquark = — ( ^ M/c Tlf c + /i e n e J . (18) 



fc 



For the NJL model we use two sets of parameters. They have been derived from those used in refs. [15| and [16| by 
neglecting the 't Hooft interaction. To do that we have followed the procedure proposed in [l7|. Namely, keeping A 
and m fixed we have varied G and m s in order to obtain M = 367.6 MeV and M s = 549.5 MeV at zero temperature 
and density. The resulting parameter sets are given in Table |TTJ 

In order to derive a quark matter EOS from the above formulae it is necessary to impose a suitable number of 
conditions on the variables {/i/ c }, /j, e , cr, ct s and A. Three of these conditions are consequences from the fact that the 
thermodynamically consistent solutions correspond to the stationary points of il with respect to a, a s , and A. Thus, 
we have 

dn/da = , dn/da s = , dn/d\A\ = 0. (19) 

To obtain the remaining conditions one must specify the physical situation in which one is interested in. In many 
astrophysical applications considered in the literature quark matter in /3-equilibrium was analyzed. In such a case 
chemical equilibrium is maintained by weak interactions among quarks, e.g. d ^> u + e~ +v e , s «-> u + e~ +v & , u+d «-> 
u + s. Moreover, for neutron stars older than a few minutes, neutrinos can leave the system. Thus, lepton number is 
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not conserved and we have four independent conserved charges, namely the electric charge uq = |n u — ^n d — ^n s — n e 
and the three color charges n u , n d and n s . Instead of n u , n d and n s , the linear combinations n = n r + n g + 715, 
713 = n r — n g and ns = -^=(n r + n g — 2nb) are often used. Here n is the total quark number density (i.e. n — 3ns) 

and t%3 and ns describe color asymmetries. The four conserved charges {rij} = {n, 713, ns, uq} are related to four 
independent chemical potentials {fJ,j} = M3, «s, mq} such that rij = —dQ/dfij. The individual quark chemical 
potentials (j,f c are given by 



fife = U + VQ 



+ /i3(A3) cc + U8(A 8 ) cc (20) 



where, as before, and A, are the Gell-Mann matrices in flavor and color space respectively. The electron chemical 
potential is fi e = — (Xq, thus we have 

«dc = fJ-sc = Hue + Me ( 21 ) 

for all colors c, which are the /3-equilibrium conditions. For electrically and color neutral matter we have also the 
conditions: 

dfl dfl dfl 

n Q = --— = , n 3 = - — = , m = -j— = . 22 
ohq on 3 d/i s 

Employing the above conditions, the system can be characterized by two independent variables, e.g. (T, ub) or (T, P). 

The conditions given in Eqs. (|21|22p are extensively employed to describe quark matter in /3-equilibrium. However, 
in this paper we deal with Just deconfined quark matter which is temporarily out of /3-cquilibrium. As already 
emphasized in [l|, 0, 0, 0, 0> H 0] the appropriate condition in this case is flavor conservation between hadronic and 
deconfined quark matter. This can be written as 

Yf = Yf f=u,d,a,L (23) 

being Yj 1 = nj /n^ and Y® = nj /n§ the abundances of each particle in the hadron and quark phase respectively. 
In other words, the just deconfined quark phase must have the same "flavor" composition than the /3-stable hadronic 
phase from which it has been originated. Notice that, since the hadronic phase is assumed to be electrically neutral, 
flavor conservation ensures automatically the charge neutrality of the just deconfined quark phase. 
The conditions given in Eq. (|23|) can be re-written as follows 

n d = £ nu, (24) 
n s = T) n u , (25) 

3n e = 2n u - n d - n s , (26) 

where u, is the particle number density of the i-species in the quark phase. The quantities £ = Yj 1 /Y^ and 
r] = Yj 1 JY^ are functions of the pressure and temperature, and they characterize the composition of the hadronic 
phase. These expressions are valid for any hadronic EOS. For hadronic matter containing n, p, A, S + , S , E - , S~ 
and 3 , we have 

^ _ n p + 2n n + n A + n s o + 2n s - + n E - 

2n p + n n + tia + 2rt s + + n s o + n s o ' 

riA + n s + + n s o + n E - + 2n H o + 2n s ~ 

rj = . (28) 

2n p + n n + ha + 2n s + + n s o + n s o 

Additionally, the deconfined phase must be locally colorless; thus it must be composed by an equal number of red, 
green and blue quarks 

n r — n g — rii,. (29) 
Also, ur, ug, dr, and dg pairing will happen provided that |A| is nonzero, leading to 

n ur = n d g, n U g = n dT . (30) 
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In order to have all Fermi levels at the same value, we consider [6| 

Tt U g — n ur (31) 
n s b = n sr . (32) 

These two equations, together with Eqs. (|29|) and ()30j) imply that n ur = n ug = n^r = nd g and n sr = n sg = n s i,. 

Finally, including the conditions given in Eq. (|19j) we have 12 equations involving the 13 unknowns (a, a s , |A|, fi e 
and {/i/ c })- For given value of one of the chemical potentials (e.g. /U ur ), the set of equations can be solved once the 
values of the parameters £, rj and the temperature T are given. Instead of \i ur , we can provide a value of the Gibbs 
free energy per baryon g quar k or the pressure P and solve simultaneously Eqs. (|2"4")) - (j52"|) together with Eq. (Tl§j) in 
order to obtain a, a s , |A|, /j, e and {/i/ c }- 

The above conditions represent a state that fulfills all the physical requirements of the just deconfined phase, i.e. 
it is color and electrically neutral, and it has a "flavor content" determined by the parameters £ and r\ (both related 
to the hadronic phase through Eqs. ([27]) and (l28|)). 

The set of twelve equations can be summarized as follows Q 

_ 2-£ -o^l _ 2r > - 

flub — * , /. T^ur i Tldb — ^ ^ , /- *^ur-, Tlsr — ^sg — Tlsb 



^ _|_ £ 5 1 ~l~ £ 1 ~f~ £ 



together with equation (fI5|> . 



2(2- 



(33) 



IV. DECONFINEMENT TRANSITION IN NEUTRON STAR MATTER 



For simplicity, the analysis that follows will be made in bulk, i.e. without taking into account the energy cost due 
to finite size effects in creating a drop of deconfined quark matter in the hadronic environment. Finite size effects on 
the nucleation of color superconducting quark bubbles in a cold deleptonized hadronic medium have been recently 
analyzed employing the MIT bag model in the description of the quark phase [2|. As a consequence of the surface 
effects, it is necessary to have an overpressure with respect to the bulk transition point. However, since this effect is 
not very large we leave its consideration for a future more detailed work. 

In order to determine the transition conditions, we apply the Gibbs criteria, i.e. we assume that deconfinement will 
occur when the pressure and Gibbs energy per baryon are the same for both hadronic matter and quark matter at a 
given common temperature. Thus, we have 

9h =g q , Ph = P q , T h = T q , (34) 

where the index h refers to hadron matter and the index q to quark matter. The results are displayed in Fig. [T] where 
we show the Gibbs free energy per baryon g at the transition point as a function of the temperature. The hadronic 
phase is described by the three parameterizations of the EOS given in Sec. II (see Table I). For the quark phase we 
adopt the NJL model described above (using flavor conservation conditions) with the two parameterizations given in 
Table II. If hadronic matter has a temperature Th and a Gibbs free energy per baryon lying to the left of a given 
curve, then the deconfinement transition is not possible for that parameterizations of the EOSs. In the right side 
region of a given curve the preferred phase is deconfined quark matter. 

Let us now examine the consequences of the above results on the structure of compact stars. Stars containing quark 
phases fall into two main classes: hybrid stars (where quark matter is restricted to the core) and strange stars (made 
up completely by quark matter). It is expected that both kinds of stars cannot exist simultaneously in Nature, but it 
is not know which one would be realized (if any). This depends on whether the energy per baryon of /?-equilibrated 
quark matter at zero pressure and zero temperature is less than the neutron mass (the so called "absolute stability" 
condition [3). Analysis made within the MIT bag model shows that there is a room in the parameter space for the 
existence of strange stars. Moreover, pairing enlarges substantially the re gion of the parameter space where /3-stable 
quark matter has an energy per baryon smaller than the neutron mass jl9l . [20| . Although the gap effect does not 
dominate the energetics, being of the order (A//x) 2 ~ a few percent, the effect is substantially large near the zero- 
pressure point (which determines the stability and also the properties of the outer layers and surface of the star) . As 
a consequence, a "CFL strange matter" is allowed for the same parameters that would otherwise produce unbound 
strange matter without pairing [19(. However, within the NJL model, the strange matter hypothesis is not favored, 
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FIG. 1: For each temperature we show the Gibbs free energy density per barion g of the deconfinement phase transition. The 
hadronic phase is described by the GM1, GM4 and GM5 parametrization of the EOS given in Sec. II (see Table I). For the 
quark phase we adopt the two parameterizations of the NJL model given in Table II. If hadronic matter has a temperature Th 
and a Gibbs free energy density per barion g^ lying to the left of a given curve, then the deconfinement transition is not possible 
for that parameterizations of the EOSs. If the point (gh,Th) lies to the right of the curve the preferred phase is deconfined 
quark matter. 



at least for the most accepted parameterizations of the EOS [17]] • Thus, stars containing quark phases are believed 
to be hybrid stars within the NJL model. 

At a given temperature there is a univocal relation between the mass of a compact star and its central pressure (or 
equivalently the Gibbs free energy per barion at the center of the star). Thus, we can employ the results given in Fig. 
[T]to calculate the critical compact star mass M cr above which they should contain a quark core. With this purpose, 
we integrate the Tolman-Oppenheimer-Volkoff equations of relativistic stellar structure employing the hadronic EOSs 
given in Tabic I, and identify the mass of the pure hadronic star for which the Gibbs free energy per barion g^ at 
the center is equal to the critical g given in Fig. [1] This is called critical mass M cr because pure hadronic stars with 
Mh > M cr are very unlikely to be observed, while pure hadronic stars with Mh < M cr are safe with respect to a sudden 
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FIG. 2: Mass-radius relationship for hadronic stars at zero temperature with the GM1, GM4 and GM5 EOS. For a given 
parametrization of the NJL model, stars with a mass above the corresponding point given in the figure are hybrid stars 
containing quark cores. The mass-radius relationship for these hybrid stars is not shown in this figure. 



transition to quark matter. The results are shown in Fig. [2] for neutron stars at zero temperature where we show 
the mass-radius relationship for hadronic stars and indicate the critical mass for the two selected parameterizations 
of the NJL model. 

In the first seconds after their formation in a core collapse supernova explosion, neutron stars may have temperatures 
up to ~ 50 MeV and a large amount of trapped neutrinos. Thus, the results presented in Fig. [I] are not appropriate 
for the analysis of the critical mass in proto-neutron star conditions since they do not include the effect of trapped 
neutrinos. However, trapped neutrinos increase considerably the critical density p cr for the transition to quark matter 
[|[. Thus, it is possible that the transition is strongly inhibited in the initial moments of the evolution of neutron 
stars About one minute after its birth, there are almost no more trapped neutrinos in neutron star matter, while 
the temperature is still high (up to ~ 10 MeV [HI). Thus, the results presented in Fig. Q]are valid for neutron stars 
older than ~ 1 minute. On the other hand, as apparent from Fig. [TJ there is little variation in the critical Gibbs free 
energy per barion for temperatures below ~ 10 MeV. Thus, the critical masses presented in Fig. [2] are also valid for 
neutron stars older than ~ 1 minute (in spite of being calculated considering hadronic matter at zero temperature). 
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V. DISCUSSION AND CONCLUSIONS 

In this paper we have analyzed the deconfmement transition from hadronic matter to quark matter and investigated 
the role of color superconductivity within the Nambu-Jona-Lasinio model. The study presented here is relevant for 
neutron stars older than ~ 1 minute, when there are almost no more trapped neutrinos in neutron star matter. 

For the hadronic phase we have used three different parameterizations of a non-linear Walecka model which includes 
the whole baryon octet and electrons. One of them is the standard paramctrization GM1 given by Glendenning- 
Moszkowski [9]. This parametrization employs "low" values for the relative meson- hyperon couplings x a , and x p . 
Larger values of these couplings make the EOS stiffer and, as shown in Table I, increase the value of the maximum mass 
of hadronic stars to values above 2 Mq. These large values may be relevant in connection with recent measurements 
of highly-massive neutron stars which give one of the most stringent test on the overall stiffness of dense matter EOS. 
However, such measurements still have to be taken with caution (see discussion given in [HI)- For example, the 
mass of the pulsar J0751+1807 was corrected from M = (2.1 ± O.1)M [H| to (1.26 ± O.14)M as new data became 
available [24| . There is also a series of measurements of extremely massive pulsars in globular clusters, where just 
the periastron advance has been determined but not the inclination angle of the orbit (2f| [2(| [27| • For the pulsar 
PSR J1748-2021B a mass of (2.74 ± O.21)M is reported by using a statistical analysis for the inclination angle 26]. 
Also, recent measurements of post-Keplerian orbital parameters in relativistic binary systems containing millisecond 
pulsars give evidence for the existence of highly-massive compact stars. For example, the compact star associated to 
the millisecond pulsar PSR B1516+02B in the Globular Cluster NGC 5904 (M5) has a mass M = 1.94±g;^M© (1 a) 
[27] ]. Other constraints for the mass and radius have been obtained studying redshifted spectral lines extracted in 
the aftermath of an X-ray burst of the Low Mass X-ray Binary EXO 0748-676 in 2000 [28[ . A model analysis of the 
X-ray burst led to rather tight constraints for the mass and radius of the compact star of M > (2.10 ± O.28)M and 
R > (13.8 ± 1.8) km [2!| which are based on the redshift measurement of (28[. However, a detailed multiwavelength 
analysis concluded that the mass of the compact star is more compatible with 1.35M than with 2.1M [30]. Moreover, 
follow-up observations of another burst in 2003 [3l[ could not confirm the spectral features seen in the burst spectra 
presented in [28|. While most of these measurements would need further confirmation, it is worth also exploring 
hadronic models that can produce stellar configurations with masses above 2 M (see the parameterizations GM4 
and GM5 in Table I). 

Employing the results of Fig. 1 and integrating the Tolman-Oppenheimer-Volkoff equations of relativistic stellar 
structure for the hadronic EOSs given in Table I we have calculated the critical compact star mass M cr above which 
they should contain a quark core. Pure hadronic stars with Mh > M cr are very unlikely to exist, while pure hadronic 
stars with Mh < M cr are safe with respect to a sudden transition to quark matter. Notice that the critical mass is 
defined here in a different way as in Refs. [H [7L l32l]. In the here- presented bulk analysis, the transition begins when 
the Gibbs conditions AP = Ph — P q = and Ag = g^ — g q = are verified at a given temperature. However, for the 
nucleation of finite size bubbles, it is necessary to have an overpressure AP > with respect to the bulk transition 
point due to surface and curvature effects. For a given overpressure there is a probability (and a corresponding 
nucleation time) to nucleate a quark bubble due to quantum or thermal fluctuations. In 0, 0] the critical mass 
was defined as the value of the gravitational mass of a hadronic star for which the nucleation time due to quantum 
fluctuations is equal to one year. Calculations of the critical mass done within the frame of the MIT bag model show 
that surface effects are strong for values of the Bag constant B smaller than ~ 100 MeV fm -3 That is, within 
the MIT bag model, the effect of the surface tension a is strong for strange stars but it is small if the parameters of 
the EOS correspond to hybrid stars. For the parameterizations of the NJL model employed in this work only hybrid 
stars are allowed. Thus, finite size effects are not expected to introduce qualitative modifications in our results. 

Employing set 1 for the Nambu-Jona-Lasinio model (see Table II) we find that there are not large variations in the 
critical mass for the three different parameterizations of the hadronic matter equation of state. As seen in Fig. 2 the 
critical mass is in the range 1.65-1.80 M for set 1, which is not a large difference if we consider the larger variation 
in the maximum mass of hadronic stars within the three parameterizations (M max between 1.78 and 2.35 M®, see 
Table I). A similar result is found employing set 2 for the Nambu-Jona-Lasinio model, i.e. the critical mass is in the 
range 1.45-1.60 M (see Fig. 2). This is in coincidence with previous results obtained for hybrid stars within the 
MIT Bag model £7} ■ As mentioned above, in the case of the MIT Bag model stars containing quark phases are strange 
stars for low values of the Bag constant B, and hybrid stars for sufficiently large values of B. For the values of B 
corresponding to strange stars the critical mass may vary essentially from zero to near the maximum mass of hadronic 
stars, depending on the value of other parameters such as the superconducting gap A and the surface tension a. 
However, for the values of B corresponding to hybrid stars the critical mass is always close to (but smaller than) the 
maximum mass of hadronic stars, rather independently of the value of other parameters [7j. The parameterizations 
employed here of the Nambu-Jona-Lasinio model allow only for the existence of hybrid stars and show the same 
characteristics of the critical mass, allowing the existence of a mixed population of compact stars (pure hadronic up 
to the critical mass and hybrid above the critical mass). 
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